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Abstract. Renormalization of non-magnetic and magnetic impurities due to electron 
double occupancy prohibition is derived analytically by an improved Gutzwiller 
approximation. Non-magnetic impurities are effectively weakened by the same 
renormalization factor as that for the hopping amplitude, whereas magnetic impurities 
are strengthened by the square root of the spin-exchange renormalization factor, in 
contrast to results by the conventional Gutzwiller approximation. We demonstrate it 
by showing that transition matrix elements of number operators between assumed 
excited states and between an assumed ground state and excited states are 
renormalized differently than diagonal matrix elements. Deviation from such simple 
renormalization with a factor is also discussed. In addition, as related calculation, we 
correct an error in treatment of renormalization of charge interaction in the literature. 
Namely, terms from the second order of the transition matrix elements are strongly 
suppressed. Since all these results do not depend on the signs of impurity potential or 
charge interaction parameter, they are valid both in attractive and repulsive cases. 
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1. Introduction 

In this paper, we discuss renormalization of impurities due to strong electron correlation. 
Such renormalization may be intuitive in the case of the Hubbard model, where each 
site has onsite electron repulsion. Namely, sites with higher potential energy have lower 
electron occupancy, and consequently have less chance of double occupancy. Then, 
the total energy loss from the impurity potential and the repulsive interaction should 
be more uniform than in the system without the electron correlation; we can call it 
renormalization of impurities. However, when the repulsion is very strong, we need to 
consider much smaller energy scales. That is, if double occupancy does not occur, the 
above argument cannot be applied, and thus renormalization of impurities within the 
lower Hubbard band is not so trivial. 

When electron double occupancy is prohibited at every site, a system with quite 
densely packed electrons has a good chance to have one electron with spin up or down 
at each site. If the system has a tendency toward phase separation, small perturbation 
by an impurity may produce a large effect to separate a system into hole-rich regions 
and electron- rich regions; it may appear in close vicinity of the half filling in the t-J- 
type models, where effective hopping is negligibly small compared to effective exchange 
interaction. In contrast, what we focus on in this paper are systems not that close to 
the half filling or systems with relatively weak exchange interaction. Then, electrons 
are more mobile. Near the half filling, since there is little freedom left to change charge 
distribution and sudden spatial change of particle number distribution around impurities 
is not favorable for the kinetic energy, non-magnetic impurity potentials may have little 
effect on low-energy eigenstates and only shift their eigenenergies quite uniformly. In 
other words, impurity potentials can be renormalized by electron correlation even within 
the lower Hubbard band. 

In previous papers [HE], such renormalization of non-magnetic impurity potentials 
was investigated numerically. That is, (i) to estimate perturbation from an impurity 
potential, the variational Monte Carlo method was applied to calculation of its 
matrix elements with respect to assumed excited states in the uniform systems; 
(ii) inhomogeneous systems with an impurity or impurities were investigated by a 
Bogoliubov-de Gennes equation with the double-occupancy prohibition treated by a 
kind of mean-field approximation called the Gutzwiller approximation (GA) generalized 
to inhomogeneous systems. 

Both of (i) and (ii) manifested strong renormalization of the impurity potential, 
and its renormalization factor (ratio between corresponding quantities in systems with 
and without the double-occupancy prohibition) seems approximately proportional to 
g* = 2x/(l + x), which is the renormalization factor of hopping amplitude obtained by 
the GA as a function of hole concentration x. Since the double-occupancy prohibition 
inhibits hopping, g l is less than unity and goes to zero as x — > 0. To explain the impurity 
renormalization factor close to g l , we pointed out the similarity between the impurity 
potential and the hopping in the real space, i.e., the Fourier-transformed impurity 
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potential has the form of hopping in the A;-space. If electrons are densely packed in 
the lattice, it must be difficult even in the fc-space to hop from A; to a different k'. 

However, it is a speculation and may not be trivial because the double occupancy 
is prohibited in the real space rather than in the fc-space. In addition, we do not 
really know how general the numerical results are because the calculation was done only 
for limited parameter sets. To complement this argument, an analytic approximation 
is adopted in this paper, namely, (i) is redone using the GA to derive dominant g l 
dependence and deviation from g t explicitly In fact, however, the conventional GA 
[31 H] fails to derive this renormalization. It compares mean weights of configurations 
relevant to operators of interest with and without the electron repulsion in calculating 
the renormalization factors. Then, the renormalization factor for the particle number 
operators is actually unity, i.e., they are not renormalized. The spin rotation invariant 
slave-boson mean-field theory [5j is known to be equivalent to the conventional GA; the 
saddle-point approximated boson fields play a role of the weights in the GA. Therefore, 
we speculate that it may have the same problem as the conventional GA. In addition, 
we believe that the slave-boson mean-field theory with only one boson often used for 
the t-J model can be even less accurate because it does not yield renormalization of the 
exchange interaction, which may be an artifact from the lost boson hard-core property. 

Let us recall that the GA corresponds to taking the leading order of the Wick 
expansion with respect to the intersite contractions of creation/annihilation operators 
[SI El- In fact, the weights of configurations in the conventional GA are likely to be 
calculated with the focus only on the lowest order; apparently it breaks down when 
the lowest order vanishes or when the next lowest order is of interest. An example is 
a particle number operator as shown in this paper. Although the lowest order is the 
average particle number, when we discuss transition matrix elements with excited states, 
this lowest order does not contribute, and the next lowest order is relevant. We will 
demonstrate that off-diagonal matrix elements between an assumed ground state and 
excited states as well as between different excited states are renormalized differently 
than diagonal matrix elements. 

Furthermore, by slightly modifying the non-magnetic impurity, i.e., by subtraction 
between up- and down-spin particle number operators, we also consider a simple 
magnetic impurity. In this case, the direction of the renormalization is reversed, namely, 
the impurity is strengthened by the electron correlation in contrast to the non-magnetic 
impurity. It must be physically reasonable because electron repulsion increases single 
occupancy. As calculation related to the non-magnetic impurity, renormalization of 
charge interaction is discussed to correct an error in its treatment in the literature. 
That is, terms relevant to the mean-field approximation are actually the second order 
of the transition matrix elements, and they are weakened by very small renormalization 
factor (g 1 ) 2 although treated usually as not being renormalized. 
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2. Model 



What we have in mind is i-J-type models with impurities, namely, 



H = Pr 



^ tij C ia C jo- + Jij ( Sj • Sj jiiflj 
i,j,cr i,j 



imp 



where c|, 

a, and Sj is the spin operator at site i. In addition 



, (T v Cj (T ) is the creation (annihilation) operator of the electron with site i and spin 



Hi = riff + n^, 



a — 1 ^%(j^%(7' 



(2) 



Gutzwiller projection operator P G = Yli(l —htffii±) prohibits electron double occupancy 
at each site and represents strong Coulomb repulsion. In this paper, we do not use any 
explicit form of tij and although they are implicitly included in assumed variational 
ground/excited states. Our main focus here is on the impurity term H iiap . 

In sections El HI and our target is renormalization of a single non-magnetic 5- 
function impurity potential located at i = I, 



imp 



Vihi = V/(n /t + fin) 



Then, in section [6], we discuss renormalization of a simple magnetic impurity, 



imp 



-hjSt = -y(n Jt -%) 



(3) 



(4) 



In addition, the focus in section [7] is not on H imp but on charge interaction hihj in 
Hamiltonian (JT]). 



3. Non-magnetic impurity renormalization 

Let us start from a uniform system without impurities. A basic idea of variational 
theories is that the ground state of the t- J-type models may be something similar to 
the BCS superconducting state 

i* ) =n( w *+ v *4t c -H) ( 5 ) 

k 

but somewhat modified by the electron correlation. Simple variational wave functions 
adopted by most of analytic theories have a form of PqI^o) with something to control 
the particle number. One way to control it is to use projection P N to fixed particle 
number ivjj] Another is to attach fugacity factors to the projector, namely, 

|*> = Pi* Q ), P = UPi, P^^X^il-n^). (6) 

i 

The latter is adopted in this paper. The reason to control the particle number is that 
Pq changes the average particle number of |^ ) because states with a larger particle 

X Many different |^o) correspond to under the projections. For example, exp(ATV) with N the 
total particle number operator is constant under P/v, and thus exp(AiV)|5'o} is equivalent to |^o)- 
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number have more chance to be projected out [SJ. Since the GA relates expectation 
values before and after the projection, 

(6)0 = <*o|0|* >, (6) = (7) 

for some operator O, usually it is not convenient if and |^ ) are totally different, 
e.g., if l^o) nas a more-than-half filled electron bandj§| Although our main interest here 
is perturbation from the uniform state, most of derivation in this paper is valid also for 
inhomogeneous systems, and thus we prefer to keep general expressions with site and 
spin indices throughout the paper, e.g., 

n ia = (n ia ) , Hi = (hi) = n it + n iir . (8) 

However, we use = (c\^c\^)q = (c| .c] .)o = (4t c il)o to avoid making formulas too 
lengthy. 

Although choice of the fugacity factors is not unique especially in inhomogeneous 
systems [7J, yet it is convenient to define 

\„ = (9) 

1 - m 

because it satisfies 

(nia) ~ (h ia )o, (10) 

for any i and a P, [7], neglecting terms of the "fourth order". Here, and throughout 
this paper, if not specified, "ra-th order" represents n-th order with respect to intersite 
contractions such as (q - c j'o-)o an d ( c 4 c it)o with i 7^ j. Note that (O) of any O can be in 
principle calculated by the Wick theorem, which yields many such intersite contractions. 
High order terms may be neglected by recalling that onsite contractions are larger than 
intersite contractions. The GA corresponds to taking the leading order only, e.g., 

(^ 2 )o « n<^ 2 >o, (ii) 



( p i)o = (1 - "itX 1 ~ n a) + -V n it(l _ n a) + -^4^4(1 - n n) 
_ (1 ~ n it)(i - n a) 



(12) 



l -Hi 

The terms neglected in the approximation in (fTTj) are of the fourth order because the 
second order terms cancel out when X ia is defined as (jSJ) [SI [7] . Let us show it explicitly 
with a notation to treat c 1 and c together, 

C ia = C iai C ia ~ C i° ') (13) 

§ The variational Monte Carlo method does not have such restriction. For example, local magnetic 
moments before and after the projection are different in general, and the chemical potential in a 
variational mean-field Hamiltonian is a variational parameter under Pn rather than a parameter to 
control the particle number. 
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by considering contractions between Pf and operators at some site(s) j,f ^ i, 

(Pi cT j'a ,C ja)o — (-^ 2 )o(cJ' CT 'Cj CT ) 

+ (-(^0^)0(^0^)0 + (ci t cJ <r )o<4 r c^) ) [(1 - ni|.) - Ai t (l - n 4 ) + A4TI4] 

+ (-(c4 c J'^')o( c 4 c Po + ( c 4 c ^)o( c 4 c I' CT ')o) K 1 - n it) + A *t n *t - M 1 ~ n *t)] ( 14 ) 
for arbitrary t, r f , a and a'. Then, the quantities in the square brackets vanish. 

We assume that |\&) is a good variational ground state, and that the excited states 
are well represented by projected quasiparticles 

|A»> = P ^ o) * PlU * o) (15) 

where jks are quasiparticles for \^f?o), namely, 

7fcf = u t c U ~ v *k c -kl, 7-J4 = v k c lt + u kC-ki- (16) 
For the denominator of \ks), we have used approximation (^o|7fc<T-P 2 7L-|^o) ~ (-P 2 )o 
[3 [9] , and errors from this approximation are of the second order. 

By switching on the impurity potential, these excited states should be mixed by 
matrix elements 

^ s Ws) „ ^^ L K (17) 

with ]Vl the number of sites. The limit of the half filling can be exactly evaluated; 
A —7- 00, PhiP —7- PP, and thus Vk',k/N\, —> (k's\ks) = 5k>k- According to the 
BCS theory, V^ s = (jk'sniTks)^) = Uk<u* k — Vk<v* k . In the previous paper [2], the 
author noted that Vk\k is not renormalized with the conventional GA [I] because it 
originally comes from a particle number operator. However, more careful analysis here 
will show that, although the diagonal matrix elements of the particle number operators 
are not renormalized [eg., see f lTUj) ]. their off- diagonal matrix elements with respect to 
the projected quasiparticle excited states are renormalized. 

The Wick expansion of ( r yk'sPni (7 P r yi s )o yields many terms, and some terms contain 
onsite contraction of hi a at the center as h Ia — > n Ia , and the others do not. Let us 
separate these two groups of terms, 

(-fk'sPn Ia P-fl s }o = n Ia (-f k 'sP 2 lL}o + (lk' s P 2 (nia - ™/ahL)o- (18) 

The first term is proportional to (k' "f \k and vanishes when k 7^ k' . Namely, we can 
only consider the second term. 

Let us first take only in the impurity potential term. Since the GA is carried 
out in the real space, the k representation should be inverse Fourier transformed into 
the real space representation. Namely, what we should calculate is {c^i a iPhi^Pc T i( ^Q. Let 
us first take the case of % ^ I, i' ^ I and i 7^ i', which makes dominant contribution to 
Vk'k- After using Pihj^Pi = A/t-n/f (1 — n/4.), we take onsite contractions for all the sites 
except i, i' and I of the numerator neglecting fourth-order terms, 
(cl! a ,Pn n Pcl) _ hM>«>Pi> n^{l-n n ) P^ a ) 



(P 2 )o ' (P§) (P?)o(P?)o 



(19) 
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Then, let us work on sites i and i', 

Pfc\ a = X ia (l - hit) C j ia , PiC ia = [(1 - n id ) + \iahia\Cia. (20) 

For the moment, we take the onsite contractions for ia and i'a' neglecting intersite 
contractions between / f or / 4- an d them; the terms neglected here are of the third 
order and will be calculated in the next section. Accordingly, using 

X i(7 (l - His) _ (1 - n i9 ) + XiaUiv _ 

(P?)o ~ {Pf)o ' { } 

f|T9|) can be approximated as 

(cJ,' a ,Ph It PcJ a ) _ A Jt (cT/ g/ w Jt (l - h n )c T ia )o 

(P 2 )o ~ (Pf)o ■ { } 

It is convenient to define mean-value-subtracted operators here, 

TT'ia = ^ia (23) 

Consequently, we obtain 

~ (c i , (T ,n n c i(7 )o - — (c^n^Oo- (24) 



(p2)o ™^,™u 

Here, the first term and the second term in the r.h.s. are from the onsite contraction of 
1 — nil an d n ih respectively; from the residual operators (n^ and 1 — hi±, respectively), 
their mean values are subtracted to cancel their onsite contraction. 

For the moment, we neglect deviation from (1241) for any i and i', which will be 
discussed in the next section. Then, it is straightforward to Fourier transform back, 

7^7 ~ (7fc'sWL>o - z (ik'sniaUo- (25) 

The formula for fij^ is obtained by exchanging "j" and J, at site /, and these formula 
represent that h Ia is renormalized into n Ia — n l5 n Ia /{\ — n l5 ). 

In fact, the derivation above is valid also for inhomogeneous systems by replacing 
7fc s with Bogoliubov quasiparticles j£. A difference is that the orthogonality of the 
Gutzwiller-projected Bogoliubov quasiparticle states is only approximately satisfied 
[7], i.e., errors from the GA can be larger than those in uniform systems. The 
renormalization of fij in inhomogeneous systems is obtained by summing up and 
n n for i + £', 

(jePfuP^o I r t _ t _ \ t \ / 0fi \ 

T u , r ~ \1t> + J Tt) > (26) 



where 



9l - (27) 
1 - ni 



'LIT 



is the Gutzwiller renormalization factor for the hopping amplitude. 

Returning to our main target, i.e., the non-magnetic uniform system, we can set 

9is = 9isi then 

Vv, k = (k'slnjlks) « g\ t {u w ul - v k ,v* k ) = gj s V™ s , (28) 
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which is exactly the same as the speculation in the previous paper [2] consistent with the 
numerical results for several /c-points by the variational Monte Carlo method, i.e., the 
renormalization factor is close to g l and insensitive to model parameters. The important 
point here may be g l appears only after summation of up and down spins, + n^, 
which is a difference from the hopping amplitude renormalization in the real space. 

According to the conventional GA [4J, what is renormalized is an operator rather 
than its matrix elements, and thus diagonal and off-diagonal matrix elements have the 
same renormalization factor. In fact, however, what is renormalized should be matrix 
elements rather than operators, and diagonal and off-diagonal matrix elements with 
respect to excited states can have different renormalization factors as demonstrated 
above. 

By exactly the same procedure as above, transition matrix elements between the 
variational ground state and projected two-quasiparticle excited states can be also 
calculated. Corresponding to (|26|) . 

i = u 1 = = ' « (Mia* fe t nj t + g^m) |*o>. (29) 
4. Corrections to the simple g l renormalization 

In the cases of i — I ^ i! , il — I ^ i and i — i' — I, we obtain formulas equivalent to 
( I24p . However, for i = %' ^ J, we have 

(c£,Pn Jt PQ _ A 7t (n Jt (l - n 7 jj cl' a ,P?c T i(T ) 

(p 2 )o ~ (PMPih ' K } 

where cJ a ,PlcJ a can be explicitly written as 

c ia Plc\ a = X ia (l - h is ){l - h ia ), c\ a PiC ia = [(1 - h id ) + \iania]h ia , (31) 

because the other combinations of cj^ and cj a yield zero or very small quantities. Then, 
although the onsite contraction of io with intersite contractions between io~ and / yields 
a formula equivalent to (I24p . the onsite contraction of io with intersite contractions 
between io and I additionally yields the same order of contribution. To compactly 
write them, let us define 

1 n in 



iCT (l-n lt )(l-n 4 ) 

as well as 



(32) 



Vi'a',ia = (fli'a'fl>ia)o ~ f^Va'^iai (33) 

which extracts only intersite contractions in (nit a th ia ) . Then, (fii> a >(l— rij (J ))o — niv(l — 
= -r]i'a',ia, and ((l-7W)(l-n i<T )} - (l-%/ (T /)(l-n i(T ) = Tfcv^. More explicitly, 

2 2 

Vi'a,ia = ~ (4 a c ia ) , r]i^ ia = (4<x c l)o • (34) 
Using these notations, 

MWO. _ _ / _ ^ \ + Kl{cLcl)o I _ Jtt_ ) . (35) 



/ D2\ ' 1 '/-I^.MT I I 2CT \ 1(7 1(7 1 U I 7J|,W 
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By summing up and n^, 

(cLP(raJt ( p 2 ^ )jPcL) ° « rbAw,* + + C(44)o (^/ t ,* + -(36) 

Since i ^ i' occurs more often than i = i', the third-order terms neglected in the 
previous section for the case of i ^ I, i' ^ I, and i ^ i' may have larger contribution 
than the newly derived terms above. Such terms are derived by taking into account 
intersite contraction including ia and i'a' . However, if intersite contractions are taken 
between ia and i'a' and the onsite contractions are taken for J, then such terms do not 
contribute as explained around ffTB"]) . Using the notation above, (120]) is rewritten as 

i>h ' i " (37) 



(^ 2 )o 

Then, for r' = —raa' (y, 4- and +1, —1 are used interchangeably) , 

V )0 

+(c T ii' c L)o [< CT (g^vm* + g^va,™) + >w (#Vfrw + ti/w 

(^ t (c} t c[,'-,)o(c/ t c[,' T ,) - ^(c} ; c[,' T ,) (c n c[,'-,)o) . (38) 

Although all the new terms in (1361) and ( 138]) contain the factors, they are 
not so simple as and inhibit the straightforward analytical transform back to k- 
representation. In other words, they cause fc-dependence of the renormalization. Since 
the ratio between the leading order and the corrections calculated in (138]) is only of the 
first order, the influence from the corrections may be larger than those in the GA for 
the real-space hopping amplitude, where the ratio is of the second order. 

Other corrections are the terms neglected in ffT5]) . We expect that they only slightly 
change the magnitude of the leading order, and that their contribution is probably not 
very important. 



5. General estimation of higher-order terms 

Let us estimate the other higher-order terms neglected above. The terms appearing 
in the Wick expansion can be classified into three groups by how to take contractions 
of hi a (l — nig): (i) Onsite contractions are taken both for / 7 and / \,. These terms 
do not contribute to Vyk as explained around (I18j) . (ii) If intersite contractions are 
taken for Ia and the onsite contraction is taken for Ia, then A/ cr n/ cr (l — hja)/ {Pj)o is 
reduced to hi a . Doing the same for \i^hi^{l — ni a )/{Pf)o yields — n/ (J n/ 5 -/(l — ni c ). 
Then their summation is g\Jhi<j. Namely, all of these terms are proportional to g\ a . 
(iii) For the other terms, intersite contractions are taken both for I J and / \,. Naive 
evaluation of these terms does not yield any explicit factor vanishing at half filling, and 
we expect that many terms cancel out each other in some way. Instead, to derive explicit 
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renormalization, let us consider such contractions for (k's\(l — n^)(l — hj^)\ks), which 
is equivalent to Vy^ for k ^ k'. Then we can replace as 
(1 -n n )(l -h n ) 1-m 

(7^ => (l-n, t )(l-n„)"'^ < 39 > 

i.e., all such terms contain ^/(l — n^) explicitly. These considerations in (i), (ii) and 
(iii) above demonstrate that Vk\k contains overall factor g\ a . 



6. Magnetic impurity renormalization 

Let us consider a simple magnetic impurity (J4|), i.e., local magnetic field is applied at 
site /. Its renormalization can be easily calculated by subtraction instead of summation 
of renormalized and hj± using formulas above. Corresponding to f[2T)j) and 

2 



\](lejP 2 1z')o(p{iP 2 le)o 



l-n^ + n n _ f 1 - n/4. + n Jt . t 
(7^mt7i>o : \lvnim)o 



(40) 



1 - Tilt ' ' 1 - »/. 

V=Z-(^S Z A) (n t = nj, (41) 



n/ 

The renormalization factor for nj-f = is (1 — n/o-)^ 1 , which is the square root of 
the Gutzwiller renormalization factor for the exchange interaction. Namely, in contrast 
to the non-magnetic impurity, the magnetic impurity is strengthened by the strong 
electron correlation. It also makes a good contrast with the unrenormalized diagonal 
matrix element (Sf) = (Sf)o (to derive this, the limit of A/f — Xii — > should be taken 
at the end starting from Aj-j- ^ Ajj,). 

In fact, also for magnetic systems (n^ ^ n ii), the factors appearing in ( 1401) are 
equivalent to those in the renormalization of the exchange interaction derived in [7J, i.e., 

(Sf sf) « (s f )„<s> + j Y. **• ^ '^f' 1 ' ) p 1 : 2 ^' ) ■ (42) 

4 <T,cr' V 1 "-to" / \ 1 / 

Although it is not explicitly noted in [7] , in this renormalization of the spin interaction, 
the first term is from onsite contractions and not renormalized (from diagonal 
matrix elements of the spin-z operators), whereas the second term including intersite 
contractions is enhanced by the renormalization factor (from the second order of the 
transition matrix elements of the spin- z operators). In fact, as shown in the next section, 
charge interaction is also renormalized in a similar manner although the direction of 
renormalization is opposite. 



7. Charge interaction renormalization 

The conventional GA [4] relates (O) to (C% for an operator O using a renormalization 
factor. By following this procedure, the renormalization factor is unity for the charge 
interaction, namely, 

? 

{flihj) « (hihjjo = UiHj + Via,ja> (43) 
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However, this approximation is correct only for the leading term n^nj and the 
renormalization factor is likely to be derived by taking only the lowest order into account. 
Using a procedure similar to that for the non-magnetic impurity, more careful analysis 
can be carried out, i.e., 



h ia (l - n^)n jV (l - n j5 i)) 




Tbi a Tlj a i T^ia ^jo" 



lioja' h 'licrjcr' 1 _ 'liaja' -i _ Viajcr' J 

-L — "Miff 1 ~~ n ja> 1 ~~ n iv 1 ~~ n ja' J 

UiHj + 9ia9%'ViT,ja' ■ (44) 



At the half filling, any state is an eigenstate of hihj with the eigenvalue unity by 
definition because every site is occupied by one electron and there is no particle number 
fluctuation, which contradicts ( I43p but is consistent with (144 p . In fact, the second term 
of r.h.s. of (|44p is the second order of fl29l) . namely, it comes from a process in which re- 
creates two quasiparticles and hi annihilates them. 

To our knowledge, every calculation in the literature on the GA is using (1431 instead 
of (jUJ) including the calculation by the author himself, and probably this error is pointed 
out for the first time here. However, this charge interaction usually does not give very 
important contribution in i- J-type models, and this correction is likely to make only 
minor modification to numerical values. Therefore, we expect that main conclusions are 
not drastically changed by this correction. Following this correction, equations in [7] 
should be modified, namely, (3<?f- — 1) and (3^ + 1) in (14) and (15) should be replaced 
by {3g!j - 9u9jj) and (3^?- + 

9ii9jj)i respectively, and derivative of gj { should be also 

considered for (16). 



8. Conclusion 



Since the Gutzwiller approximation is formulated to (almost) conserve the particle 
number at the Gutzwiller projection, one may consider that quantities related to 
particle number operators are not renormalized. However, since the particle number 
is an expectation value with respect to an assumed ground state, the constraint of its 
conservation does not restrict transition matrix elements with excited states. Our results 
here correct description by the conventional Gutzwiller approximation in the literature, 
where such renormalization factors are calculated with a focus on diagonal matrix 
elements or lowest-order terms and regarded as unity. The results in this paper are 
general and do not depend on parameters. Namely, they are valid both for attractive and 
repulsive impurity potentials and both for attractive and repulsive charge interactions. 

The Fourier-transformed impurity potential has a form of hopping in the fc-space. 
We have derived similarities and differences between this "hopping" in the /c-space 
and in the real space under real-space electron double-occupancy prohibition. As a 
similarity, they are strongly renormalized to decrease with hole concentration x, and 
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their renormalization factor is g l = 2x/(l + x) in uniform non- magnetic systems. In 
addition, the higher order terms also contain g l . It should represent that not many 
available seats to hop are left because of the electron repulsion. A difference is, 
however, (c\ a Cj a ) of each a is renormalized in the real space, whereas renormalization 
of Y^a{ c k'cr c ka) appears only after the summation over spin a = ± in the /c-space. If this 
summation is replaced by subtraction, which corresponds to a magnetic impurity in the 
real space, then the direction of the renormalization is reversed, i.e., the renormalization 
factor is larger than unity and equivalent to the square root of that for the exchange 
interaction. As another difference, the corrections to the leading order term in the k- 
space can be larger and have more complicated expression than those in the real space. 

As related calculation, renormalization of charge interaction has been also derived. 
The leading order is rather trivial and unrenormalized, i.e., it is the product of particle 
densities at the two relevant sites. The next leading order term is the second order 
of transition matrix elements of the number operators with excited states. Since 
the transition matrix elements are renormalized by g l , these second order terms are 
renormalized by (g*) 2 , namely, strongly reduced. These terms include hopping and 
pairing amplitude and are relevant to the mean-field approximation. Similar relation 
is found also in the ^-component of the exchange interaction. Namely, the leading 
order is the product of spin-z densities at the two relevant sites. The next term is the 
second order of transition matrix elements of the spin-z operator, which is strengthened 
by the electron repulsion. At the half filling, any state is an eigenstate of hifij, with 
the eigenvalue unity. In fact, (jSJ) satisfies it even in magnetic systems, which may 
demonstrate that the choice of fugacity factors by ([9]) is reasonable. Other choices of 
fugacity factors also discussed in [7] do not seem to satisfy it in magnetic systems, and 
their use is likely to be restricted in systems with small magnetic moments. 
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